WIMP dark matter, Higgs exchange and DAMA 
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In the WIMP scenario, there is a one-to-one relation between the dark matter (DM) relic density and spin 
independent direct detection rate if both the annihilation of DM and its elastic scattering on nuclei go dominantly 
through Higgs exchange. In particular, for DM masses much smaller than the Higgs boson mass, the ratio of the 
relevant cross sections depends only on the DM mass. Assuming DM mass and direct detection rate within the 
ranges allowed by the recent DAMA collaboration results -taking account of the channelling effect on energy 
threshold and the null results of the other direct detection experiments- gives a definite range for the relic 
density. For scalar DM models, like the Higgs portal models or the inert doublet model, the relic density range 
turns out to be in agreement with WMAP. This scenario implies that the Higgs boson has a large branching ratio 
to pairs of DM particles, a prediction which might challenge its search at the LHC. 



The DAMA collaboration has recently provided evidence for 
an annual modulation of the rate of nuclear recoils in their 
detector [[0], confirming at a firmer level their previous re- 
sults |Q]. Taking into account the null results of the other 
direct Dark Matter (DM) detection experiments [3D and the 
recently discovered channelling effect on the threshold energy 
in DAMA, points towards a nuclear recoil due to dark matter- 
nucleon elastic scatterings, with a spin independent (SI) cross 
section in the range (see Q]) 

3 x 10" 41 cm 2 <cf <5x 10" 39 cm 2 (1) 

and dark matter mass in the range 

3GeV<m DM < 8GeV. (2) 

These results have been already the object of various studies in 
specific models (ESS Si- In this short letter, working in the 
WIMP framework which assumes a DM relic density deter- 
mined by the thermal freeze-out of DM annihilation, we em- 
phasize the importance of Higgs exchange diagrams for DM 
mass in the range of ©. 

If the DM candidate is light, there is a limited number of pos- 
sible (2-body) annihilation channels to SM particles, i.e. to 
u,d,s,c,b quark pairs, to lepton pairs and to photon pairs. 
Depending on the model this can be done at tree or loop 
level in various ways. Annihilations through the SM Z bo- 
son may be excluded right away for it would imply that dark 
matter contributes to the Z invisible width [ 9]. The next sim- 
plest possibility at tree level is annihilation of dark matter into 
fermions through the Brout-Englert-Higgs (Higgs for short) 
in the s-channel, as in the diagram of Fig.l.a. In this case, 
only bb, cc and XX annihilations are relevant, since all other 
SM fermions have small Yukawa couplings. From the very 
same coupling between DM and the Higgs, SI elastic scatter- 
ing is induced through a Higgs in the t-channel, Fig.l.b. For 
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FIG. 1: Higgs exchange diagrams for the DM annihilation (a) and 
scattering with a nucleon (b). 

the Inert Doublet Model (IDM), a two Higgs model exten- 
sion of the Standard Model with a scalar dark matter candi- 
date id [H El H H H d, and a fortiori for a singlet 
scalar DM candidate HIIlElliSlI^ chan- 
nels of Fig.l are the only possible non-negligible ones in the 
range of ©. In more sophisticated models, such as with the 
neutralino DM candidate of the MSSM, these channels coex- 
ist a priori with many other channels, in particular with other 
intermediate Higgs scalar particle or squarks channels 1F7LI25I1. 



Both processes in Fig.l are tightly related. The Higgs bo- 
son mass dependence is the same, m^ 4 , because in both cases 
the momentum of the Higgs is negligible compared with m^. 
Furthermore the dependence in the unknown DM-DM-/1 cou- 
pling is the same. This means that, up to uncertainties in the 
Yukawa couplings of the SM fermions and in the Higgs to nu- 
cleon coupling, the only left parameter is the mass of the dark 
matter. This dependence is however limited if we take the DM 
mass within the range ©. In other words, if Higgs mediated 
processes are dominant, the ratio of the cross sections corre- 
sponding to the processus of Fig.l is essentially fixed in any 
model and there is a one-to-one relation between annihilation 
and direct detection. This has been pointed out for the case of 
scalar singlet DM in Ref. [[HI] . The ratio may however be dif- 
ferent in different models. To see this and to inquire whether 
any model may be in agreement with observations, we con- 
sider two simple scenarios, respectively with a scalar and a 
fermionic dark matter candidate. 
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For a scalar dark matter candidate, the simplest possibility is 
to introduce one real scalar singlet S, odd under a Z2 sym- 
metry G2L EE [H S [2TU2I 111 0], like in the so-called 



Higgs portal framework |20l|2JJ,|22l|23|]. In full generality the 
four following renormalizable terms may be added to the SM 
lagrangian: 

L 3 h^SdfJS - \i? s S 2 - ^-S 4 - X L H^HS 2 (3) 

with H = (h + (h + iG )/V2) T the Higgs doublet. The mass 
of S is thus given by 



2 2,^2 
m S=I Ll S + ^ v • 



(4) 



where v = 246 GeV. In this model the sole coupling which 
allows S to annihilate into SM particles and to interact with 
nuclei is Xl. For the annihilation cross section, the elastic 
scattering cross section (normalized to one nucleon) and the 
ratio, we obtain 
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where ^ c = 3(1) for quarks (leptons), v r£? / = (5 — 4m|) 1 / 2 /m i s 
is the centre of mass relative velocity between both S and 
ju r = m^m^v / (ms + m#) is the nucleon-DM reduced mass. The 
factor / parametrizes the Higgs to nucleons coupling from the 
trace anomaly, fm^ = (N\ Y^ q m q^ c l\^) — ghNNV- From the re- 
sults quoted in Ref. (26D, we take / = 0.30 as central value, 
and vary it within the rather wide range 0.14</<0.66. As 
for the Yukawa couplings, Y{ = \/2mi/v, we consider the pole 
masses = 4.23 GeV, m c = 1.2 GeV and m T = 1.77 GeV 
(we neglect the effects of the running of the Yukawa couplings 
which are expected quite moderate). 

Within the WIMP scenario, one needs a(SS —> all ff)v re i ~ 
1 pb to have the right dark matter abundance. Therefore using 
Eq.([T]), the ratio required by data is R 10 3-4 . This has to 
be compared with the value of R obtained from Eq. ©. Fig. 2 
gives both ratios as a function of ms, requiring that the relic 
density with respect to the critical density obtained is within 
the WMAP density range 0.094 < £l DM h 2 < 0.129 S Hi, 
and that c s p r and ms are in the DAMA region allowed by other 
direct detection experiments, see Fig. 1 of Ref. [4] (taking into 
account the channelling effect). As R scales approximately as 
m\ for ms > ray, it is remarkable that both values of R coincide 
around the DAMA DM mass region. The relic abundances are 
computed using MicroMegas Il29llll53ll. 
This can also be seen in Fig. 3 which gives the regions of ms 
and Xl (for m^ = 120 GeV) consistent with WMAP and the 
same direct detection constraints. Both regions nicely overlap. 
The corresponding values of the SI elastic cross section can be 
read off from Fig.4. For the central value / = 0.30 the overlap 
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FIG. 2: Values of the ratio R calculated from Eq. Q for three values 
of / to be compared with the area of R values required to match at 
2a level both WMAP relic density and direct detection constraints 
(i.e. DAMA, channelling effect included, and upper bounds from 
CDMS, Xenon, CoGent and Cresst, see Fig. 1 of Ref. H). 
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FIG. 3: For = 120 GeV, values of ms and Xl which lead to 
the WMAP result, 0.094 < £l DM h 2 < 0.129 (solid black lines), and 
which match the direct detection constraints {i.e. Fig.l of Ref. 0]), 
for the central value / = 0.30 as well as the values / = 0.14 and 
/ = 0.66. 

region covers the ms ~ 6-8 GeV range while for 0.14 < f < 
0.66 regions overlap for ms between 3.5 GeV and 8.5 GeV. 
For / smaller than 0.20 there is no overlap region. For a fixed 
value of ms, a smaller / gives a smaller detection rate. This 
might be compensated by a larger Higgs-DM coupling Xl, but 
then at the expense of a smaller relic abundance. 

In Figs. 2-4 we used mu = 120 GeV but agreement may 
be obtained for other Higgs boson mass provided the ratio 
Xhlm\ is kept fixed. Typically the required value is Xi/m^ ^ 
10 -5 GeV -2 . To keep the result perturbative, Xl < 271, we 
need that m h < 800 GeV. 

Since all the parameters are fixed, or strongly constrained, we 
may also make some definitive predictions regarding possible 
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FIG. 4: For m h = 120 GeV, \ogaf (in pb= 10" 36 cm 2 ) as a function 
of ms and \l, versus Q^dm^ 2 (black lines). For other values of m/ t 
a S p scales as 1/m^. 

indirect detection of DM, in particular through gamma rays 
from annihilation of DM at the galactic centre. The mass of 
© puts the DM candidate in the energy range of EGRET data 
and of the forthcoming one of GLAST. Fig. 5 shows the pre- 
dicted flux of gamma rays from the galactic centre for a sam- 
ple of scalar DM with parameters which are consistent with 
DAMA and WMAP. EGRET data f47|] are also shown HI. It 
is interesting that the predicted flux is of the order of magni- 
tude of the observed flux at the lowest energies that have been 
probed by EGRET. Actually, the predicted flux may even be 
larger than the observed one for some set of parameters. We 
have however refrained from putting constraints on the model 
parameters, given the large uncertainties on the abundance of 
dark matter at the centre of the Galaxy. Here we assume a 
mundane NFW profile [45], as in 111 411 . Similar predictions 
regarding the gamma flux, for a different model compatible 
with DAMA, have been done in 

A number of remarks are now in order. First of all, we should 
keep in mind that the parameter fit to detection data depends 
on precise assumptions on the abundance and velocity of the 
distribution of dark matter in our neighbourhood (see e.g. |4|]). 
Also, the cosmic relic abundance is obtained assuming a mun- 
dane, radiation dominated expansion of the universe. These 
are conservative but well motivated options. However relax- 
ing either one or the other would give other predictions and 
would likely jeopardize the model. 

On the experimental side, it is to be seen if the intriguing re- 
sults of DAMA will survive the test of time and will be even- 
tually confirmed by other SI experiments, like CRESST (see 
i.e. the discussion in [4]). Further confirmation of the chan- 
neling effect in the regime studied by DAMA would also be 
useful. Without channeling, there is no region allowed by all 
experiments ll55ll ll56ll . Regarding SD detectors, as there are 
no (non-negligible) spin dependent (SD) interactions in the 
Higgs exchange scenario, no signal should be expected. 
On the theoretical side, we definitely wish we had a better de- 
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FIG. 5: For three example values of ms and ]U2 (in GeV) taken from 
Fig. 3, flux of gamma rays from the galactic centre from the annihila- 
tion of a scalar DM consistent with DAMA, compared with EGRET 
data. 



termination of the form factor / in the Higgs-nucleon effective 
coupling (see the discussion above). 

Also on the theoretical side, inspection of Fig. 3 and 4 reveals 
that substantial tuning is necessary to obtain agreement with 
DAMA data, as one needs a quite large SI cross section. As 
already emphasized in fl9ll (see also i22h for a singlet and 
iB for the IDMlEl) the cross section of a scalar DM can- 
didate is typically small, 10 -8 /?Z?, below the range of all the 
existing direct detection experiments. For it to be observable 
in DAMA requires a large, albeit still perturbative coupling 
\X L \ ~ 0.1 - 1 between the DM and the Higgs |58|]. Simul- 
taneously, to keep the DM mass within the range © requires 
a cancellation between both terms in Eq. © at the per cent 
level, roughly. This, in our opinion, is not unbearable. Nor is 
it suprising, given the very minimal number of parameters of 
the model. Taken at face value, it implies a close connection 
between DM and the Higgs sector, thus with the mechanism 
at the origin of electroweak symmetry breaking. 
This brings us to one immediate, striking consequence of such 
DM models. A large coupling to the Higgs leads to a large 
Higgs boson decay rate to a scalar DM pair [19]. For ex- 
ample for ms — 1 GeV, Xl = — 0.2 and for a Higgs of mass 
120 GeV we get the branching ratio BR(h -> SS) = 99.5%, 
while for m h = 200 GeV and X L = -0.55 we get BR(h -> 
SS) ~ 70%. This reduces the visible branching ratio accord- 
ingly, rendering the Higgs boson basically invisible at LHC 
for nih = 120 GeV, except possibly for many years of high 
luminosity data taking. Such a dominance of the invisible 
DM channel is a clear prediction of the framework, although a 
challenging one for experimentalists for low value of nth (see 
[49] for strategies to search an invisible Higgs at the LHC). 



We now would like to emphasize that the results discussed 
here apply for any scalar dark matter model for which annihi- 
lation and SI scattering cross sections would be dominated by 
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the diagrams of Fig. 1. A simple such an instance is the IDM, 
in which case the DM candidate is the lightest component of 
a scalar doublet, odd under a Z2 symmetry. All the results 
above hold, provided one identifies the appropriate parame- 
ters. Using the notations of 111411 for instance, we get the same 
abundances and SI cross section in the DAM A range, Figs. 2- 
4, provided one replaces ms — > m# , jus — > m and Xl^^l- 
One just has to make sure that the extra components of the 
doublets, noted Ao and H ± , are heavy enough. In practice, 
this means that we need m# + ma > mz so as not to con- 
tribute to the Z invisible decay width. This simultaneoulsy 
kills the possibility of sizable co-annihilation of Ho with Ao in 
the early universe or elastic scattering in detectors through a Z 
channel. As we consider m# in the range 0, ma > mz- Fur- 
thermore we need to preclude large radiative corrections to the 
Z and W bosons, and in particular to the p parameter (see the 
discussion in 111211 ). Within the IDM, this may be obtained nat- 
urally as there is a global custodial symmetry if rriA ~ m H + , 
in which case the contribution of the extra scalar doublet to 
p vanishes (see ll3lll and the discussion in lfl6ll ). Of course, 
the extra heavy degrees of freedom Ao and // ± might show up 
eventually in high energy collisions D2I1 . something in which 
the IDM differs from the simpler singlet scalar model. 
For the sake of comparison with the scalar DM case discussed 
above, we now briefly consider the case of singlet Dirac and 
Majorana fermionic DM candidates. In the former case, we 
consider the Lagrangian 



L 3 \|/($-mo)\|/- 



(8) 



The Higgs-DM coupling may arise from the non- 
renormalizable operator \\n\fH^H/A, so that Yy oc v/A, 
and the mass of \|/ is — mo + Yyw/ y/2 (for similar models 
see i.e. (33D). For the annihilation and direct detection cross 
sections we get 



c(\\n\f^ff)v rd = n t 



Y 2 m }v 2 rel (rn^-rn}) 3 / 2 



1671 y 2 m\ 



a(\\fN — > \\fN) 



v2 2 

2k v 2 m\ 



(9) 



R 



_ E/C>(W ff)v re i _ T.fn c m 2 f v 2 rd (m 2 - m 2 f fl 2 



Both the cross sections and the ratio have different paramet- 
ric dependence compared to the DM scalar case, Eqs. dS])-©. 
Keeping all things constant, the SI independent cross section 



is smaller by a factor of m^ M /\ 2 . Also, assuming perturbative 
values of the coupling Y^, the abundance of dark matter is way 
larger than the critical density [for parameters consistent with 
CQ) and 0] because of the extra small factor mj) M /v 2 • v 2 rel 
in the annihilation cross-section compared to the scalar DM 
case. The origin of this suppression is well-known (see e.g. 
ll34ll ). It results from the fact that a fermion-antifermion pair 
in a s-wave is a CP odd state and can not annihilate into a 
scalar, like the Higgs. Hence the annihilation is both p-wave, 
oc v 2 rd , and helicity suppressed, oc m 2 DM . Similar conclusions 
hold for the Majorana case. In these fermion cases other chan- 
nels than Higgs exchange must necessarily be present in order 
to match the WMAP and DAMA observations, as for instance 
is the case in the MSSM (see the recent discussion of 10]). 
Before concluding, let us raise that the scenario discussed here 
does not shed light on the baryon-dark matter coincidence 
problem, £Idm ~ If the dark matter mass is in the range 
0, there is roughly one dark matter particle per baryon in the 
universe, ubm^^b^ result which suggests a deeper connec- 
tion between ordinary and dark matters. A few models have 
been proposed to naturally explain this coincidence (see e.g. 
^^Wi^W^^^^^k)- They tend to be so- 
phisticated in comparison with the model discussed here and, 
to our knowledge, they also tend to have cross sections for 
scattering of dark matter with nuclei that are much below the 
sensitivity of existing and forthcoming detectors. 

In conclusion, we have discussed the possible relevance of 
Higgs exchange for WIMP in DAMA mass and SI cross sec- 
tion ranges, leading to a one-to-one relation between the relic 
density and SI direct detection rate. We have emphasized 
that for scalar dark matter this relation is in agreement with 
data, provided Higgs exchange is the dominant channel in 
both annihilation and SI cross sections. We have illustrated 
this through a singlet real scalar and the inert doublet model, 
where it is naturally realized. These models are very simple, 
give definitive predictions -including for the LHC- and might 
be falsifiable within a foreseeable future. 
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